Gauge/gravity duality and thermalization of a boost-invariant perfect fluid 
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We derive the equation for the quasi-normal modes corresponding to the scalar excitation of a 
black hole moving away in the fifth dimension. This geometry is the AdS/CFT dual of a boost- 
invariant expanding perfect fluid in A/" = 4 SUSY Yang-Mills theory at large proper-time. On the 
gauge-theory side, the dominant solution of the equation describes the decay back to equilibrium 
of a scalar excitation of the perfect fluid. Its characteristic proper-time can be interpreted as a 
thermalization time of the perfect fluid, which is a universal (and numerically small) constant in units 
of the unique scale of the problem. This may provide a new insight on the short thermalization-time 
puzzle encountered in heavy-ion collision phenomenology. A nontrivial scaling behaviour in proper- 
time is obtained which can be interpreted in terms of a slowly varying adiabatic approximation. 
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Thermalization response-time and Black Hole 
quasinormal modes 

In a recent paper m, we have shown that the 
AdS/CFT dual of an expanding relativistic thermal and 
perfect fluid in A/" = 4 SUSY Yang-Mills (SYM) theory 
can be identified as a black hole (BH) moving away in the 
fifth dimension In the above paper, a holographic 

renormalization procedure Q using Fefferman- Graham 
coordinates 0] allowed to construct the gravity duals of 
a continuous 1-parametric set of 4-d stress-energy tensors 
of the gauge theory containing among others the duals 
of the free streaming and perfect fluid cases. Interest- 
ingly, the corresponding family of geometries was shown 
to possess singularities except only for the latter case, 
which happens to be a BH moving away in the fifth di- 
mension. This correspondence was shown to be valid at 
asymptotic proper-times, independently of initial condi- 
tions provided boost-invariance is preserved. This gives 
an AdS/CFT physical criterion for the emergence of per- 
fect fluid behaviour at large proper times. 

Starting from this correspondence, it is interesting to 
study the stability properties of this BH system, since 
it can bring a physical insight on the typical relaxation 
proper-times of the relativistic thermal and perfect fluid, 
which does not appear reachable from a direct strong 
coupling computation in the J\f — 4 SYM field theory. 
As was done in the static BH case |^ in order to calcu- 
late the decay time of an excitation of the system, one 
computes the quasinormal modes (QNM) of the Einstein 
equations linearized around the background geometry. In 
particular the calculation of the QNM's corresponding to 
a scalar excitation canonically coupled to the metric in 
the gravitational dual configuration 0, S 0, S] may give 
an evaluation of the thermalization time of the dual gauge 
field-theoretic strongly-coupled system after a small de- 
viation from equilibrium. 

In the present paper, our aim is to extend the analy- 



sis from the static case to the black hole moving off in 
the fifth dimension, deriving the equation and making 
the evaluation of the corresponding QNM's in order to 
estimate the thermalization decay proper-time of the rel- 
ativistically expanding perfect fiuid in the A/" = 4 SYM 
field theory. 

The plan of our study is the following. In section 2, we 
shall (re)derive the QNM equation for the static BH using 
now the Fefferman- Graham coordinates. In section 3 we 
derive the QNM equation for the gravitational dual of the 
expanding perfect fluid and give its solutions both in the 
minimally coupled scalar case and for transverse tensor 
perturbations. In section 4, we discuss the features of 
our results and their possible physical relevance for the 
thermalization puzzle of the QCD quark-gluon plasma in 
heavy-ion collisions. We conclude and give an outlook. 



Quasinormal modes of a static Black Hole in 
Fefferman- Graham coordinates 



Quasinormal modes (QNM's) formally define the re- 
sponse of a black hole state to small perturbations, for in- 
stance due to a scalar-field excitation canonically coupled 
to the metric, with incoming (absorbing) boundary con- 
dition at the BH horizon. The main observation is that 
the resulting frequencies become complex and hence the 
perturbations are expected to die out exponentially |0. 

Quasinormal frequencies have been calculated for nu- 
merous examples of static black holes in various num- 
bers of dimensions, in particular for the static planar 
AdS black hole which is the dual geometry to A/" = 4 
SYM theory at nonzero temperature. In this section 
we shall rederive the known results using the Fefferman- 
Graham coordinates which are suitable for the extension 
to the relevant non-static case. As already mentioned 0] , 
the expanding geometry has a simple form when directly 
written in these coordinates. 
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In FefFerman- Graham coordinates the metric of a 5-d 
planar BH has the form: 
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dx^ + dz^ 



(1) 

where (t, x) are the boundary coordinates, z is the fifth 
one, and zq is the location of the horizon in the bulk. 

Quasinormal modes for a scalar perturbation of a black 
hole are obtained by solving the wave equation for a 
massless scalar field 



A(j) : 



'9 



(2) 



where g^^ is the metric tensor and g its determinant in 
the background geometry J^l, assuming purely incoming 
boundary conditions at the horizon z = zq and Dirichlet 
conditions at the boundary (see e.g. 0). 

Inserting eq. iQJ into (0) , the equation for a scalar field 
(with zero transverse momentum) has then the form: 

1 (i-S) ... . . fi 



dt^{t,z) + dA^{i-^]dA{tr^ 



A separation of variables 

0(i,z) = e^"V(2) 
leads to the ordinary differential equation 
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Note that, by a change of variable 

z/zo = tanh^/^(4z*) , 
the equation takes the canonical form 



V8 e^^^' 
sinh^/2(8z* 



-uj^(t){z*) = 



(3) 
(4) 

. 

(5) 
(6) 
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from which one can determine the quasinormal frequen- 
cies UJ. 

In fact, by a further change of variable 

i^"-';^'°':)^i-2''°"'"';'''' (8) 

1 + {z/zqY 1 + tanh(4z*) ' ^ ' 

one shows that this equation can be put in the form of 
the Heun equation 

<^"+— — r <t>' + ^ 0^0 (9) 

^ z(l - z)(2 -i)^ 4z(l - z){2 -z)^ ^ ' 

which was obtained 0] starting with the conventional 
BH metric. The QNM solution of jnj dominant at large- 
time (i.e. with smallest imaginary part) is found to be 
i^staUc/^T ^ 3.1194 - 2.74667 i 



Quasinormal modes for the boost-invariant Black 
Hole geometry 

In this section we derive the equation for a scalar field 
in the background corresponding to the asymptotic ex- 
panding perfect fluid geometry. The geometry has the 
form [i|l3: 



ds' = 



(1 



-dr^ + (1 + v^) (r^dy^ + dxf) + dz^ 



(10) 

where Xi=i^2 sue the transverse coordinates, while the 
proper-time r and rapidity y are related to the longitu- 
dinal coordinates through t = rcoshy and 0:3 = rsinhy. 
w is a scaling variable 

z 

v = —. (11) 
T3 

To this background metric let us couple a scalar field 
which depends only on the proper-time r, and on z. The 
equation ^ takes the form: 

(12) 

Since the metric (fTT)|l has a simple dependence on the 
variable v let us rewrite the above equation in terms of 
r and v. The matrix of differentials is 



-dr 



drA+d. (^(l-i)«)9.0) 







(13) 



Since the perfect fluid geometry is valid in our problem 
only at large proper-time, in the scaling limit v = const 
and r — > 00 we may consistently neglect the non-diagonal 
contribution in (|13|l . Thus, the QNM calculation pre- 
serves the specific scaling property in t/z^ = of the 
perfect fluid solution of the AdS/CFT correspondence. 
The resulting equation takes the form: 

1 {1-vy 
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5^0(r, v)+T—^ dy ^(1 - ««)9„0(r, v) 



(14) 

which has some similarity with the one for the static black 
hole, once substituting variables {t,z) {t,v), see 
The noticeable difference is the additional t~'^I'^ factor, 
which leads to a nontrivial scaling in the (proper) time 
dependence. 

Performing a separation of variables (/)(t, u) = /(t)(/)(u) 
we get two decoupled equations. An important point is 
however to notice that the equation for the r-dependence 
is no longer a plane wave as in but is determined by 
the equation 



dl!{T)^-U?T-^j{T) , 



(15) 



whose solutions are linear combinations of the Bessel 
functions 
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(16) 



3 



In the large r behaviour, the region of vahdity of the 
perfect fluid geometry being asymptotic, the relevant be- 
haviour of the Bessel functions is 

/(r) ~ ei--^ , (17) 

to be compared with the plane waves (QJ of the static 
case. We will comment on the physical interpretation of 
(|17|l in the next section. 

The resulting r-independent factorized equation reads 

dv f - v'')d.^{v)) +uj^l ^\-^'}% {v) = (18) 

which, strikingly enough, is formally the same ordinary 
differential equation for the r-independent part as eq. 
for the static black hole but with v playing the role of z. 
Hence we get the same profile functions (j) of the quasi- 
normal modes and their eigenfrequencies uj. However the 
variables in the perfect-fluid case v and r are different 
from the z and t variables relevant for the static black 
hole. Indeed, restoring the z and t dependence of the 
evolving solution 4>{v)f{T) gives rise to a rather intricate 
spacetime dependence. 

It is also interesting to derive the quasinormal modes 
for perturbations of the metric which on the gauge 
theory side correspond to perturbations of the energy- 
momentum tensor. Of particular interest are perturba- 
tions of the component T^^x^- Physically their exponen- 
tial decay can be interpreted as transverse isotropiza- 
tion of the asymptotic hydrodynamic expansion of the 
plasma. On the gravity side they are interesting since, 
as shown in |lfll | for general static cases, their equation 
of motion is identical to the scalar one. In particular 
let us consider the metric component gxiX2 ^-^d form the 
quantity 

gll ^ z\-^(^^^-^gx,x, ; e^^"'^) ^ 1 ^ (19) 

where c{v, z) is the solution 1] found for the transverse 
component of the boost-invariant perfect fluid metric (see 

eq. cni)). 

We have shown that in the case of the moving black 
hole, in the large r limit, the quantity g%l^ also satisfies 
also the scalar equation of motion. 

(20) 

where the proper-time dependence has been separated 
out in the same manner as for the scalar. Consequently 
this mode of the metric has the same set of quasinormal 
frequencies. 

Conclusion, discussion and outlook 

Let us summarize our results. We consider the problem 
of computing the response proper-time of an expanding 



relativistic and boost-invariant M = A SYM perfect fluid, 
after a scalar excitation off thermodynamic equilibrium. 
In the dual gravitational geometry, which corresponds to 
a black-hole moving away in the fifth dimension, the char- 
acteristic proper-time is related to the frequency uJc with 
the smallest imaginary part among quasi-normal modes. 

We derived the corresponding equations in the scalar 
and the transverse tensor channel. The resulting equa- 
tions factorize into two decoupled equations defining re- 
spectively the explicit proper-time and scaling variable 
dependence of the modes. As in the static case, we find 
that the quasinormal frequencies have imaginary parts 
which correspond to exponential decay of perturbations 
towards the (expanding) equilibrium state. While for- 
mally the equation defining the frequencies has the same 
functional form as in the static case, they correspond to 
quite different variables namely the proper time r and 
the scaling variable v ~ z/t^/^ . 

The profile of the solution is given as a function of the 
scaling variable w, </)(«), which obeys the same equations 
as 4>{z) in the static case. By contrast with the plane 
waves of the static case however, the frequencies are re- 
lated to Bessel functions ^/tJj^z (I'-^'''^) ~^ '''^ e^*'^'^^, 
hence leading to a non-trivial scaling in proper-time. 

These results can be interpreted as follows. The tem- 
perature for the expanding perfect fluid is known to be- 
have like T ^ r~^/^ which is consistent through 
the AdS/CFT correspondence with the evolving black 
hole solution of P|. In order to understand the nontriv- 
ial scaling in proper time (|17() . we note that the static 
quasinormal freqencies are proportional to the tempera- 
ture i.e. LOstatic{T) = a ■ T, where a is a constant. In 
the expanding case we can consider an adiabatic approx- 
imation assuming that locally the temperature is fixed. 
Hence a plane wave dependence 

2. 

will give the scaling in t^/'^ seen in the overall solution. 
This adiabatic approximation deserves further study. 

A similar discussion is valid for transverse perturba- 
tions of the metric. It would be interesting to extend the 
analysis to generic metric perturbations. 

On a theoretical ground, a comment is in order about 
the well-known calculation of the viscosity to entropy 
ratio [l]| in the static BH configuration. Quasinor- 
mal modes and viscosity calculation correspond both to 
poles of specific retarded Green functions for a scalar 
(for QNM's) and metric (for viscosity) deformation us- 
ing Fourier transforms in the time variable t. This is the 
case for the static case, due to the plane- wave solution 

The difference between the scalar QNM's and viscosity 
(shear channel) is that the corresponding poles do not go 
to zero at small transverse momentum contrary to the 
viscosity case 0. This is expected, since viscosity is a 
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hydrodynamic excitation surviving at large time scales, 
while normal QNM's lead to an exponential fall-ofF, at 
least within an AdS space. This is the origin of the finite 
"thermalization response-time" obtained for the scalar 
equation. 

Now turning to the expanding geometry in proper- 
time, viscosity remains an interesting issue for future 
work, since it is not clear whether one can use the method 
of Fourier-transformed retarded Green functions to eval- 
uate it, since the Fourier transforms seem not to be di- 
rectly relevant for the proper-time dependence. An adi- 
abatic approximation could help with this problem. We 
postpone this analysis for future work. 

In a phenomenological perspective, it is instructive to 
make a parallel between our calculation in the frame- 
work of the A/" = 4 SYM fluid, and the problems dis- 
cussed about the thermalization proper-time of a QCD 
plasma formed in a heavy-ion collision at high energy 
(see ^3 and references therein). Indeed, there are some 
indication that this typical proper-time is rather short, 
which is difficult to explain in terms of initial conditions 
dominated by a weak-coupling state, and thus with high 
viscosity and a-priori long thermalization time. 

In m, we have shown that the gauge/gravity corre- 
spondence was selecting the perfect fluid solution at large 
proper-times. Since it seems that the model of a QCD 
perfect fluid could be favored by the phenomenological 
analyses, it is perhaps not completely unrealistic to look 
for some physical lessons of our present results. The main 
point is the rather strong stability of the asymptotic so- 
lution, since the response-time to a scalar excitation is 
short in terms of proper time. Combining eq. (|17|l with 
the numerical value of the dominant QNM leads to a 
proper-time damping of the form 

exp -2.7466 -ri^ (22) 

in the units where the horizon is fixed at zq = t~^/^. 

If these estimates would be also approximately valid for 
QCD, at least in the hypothesis of a deconfined phase for 
which the supersymmetric degrees of freedom would not 
play a major role, this can lead to a new point-of-view on 
the thermalization problem. Indeed, the stable, strongly 
interacting state represented by the perfect gauge fiuid, 
could act as an "attractor" during the proper-time evo- 
lution of a collision with QCD plasma formation, such 
as for heavy-ion collisions. This would give a typical 
non-perturbative mechanism, a priori complementary to 
those related to the initial perturbative conditions. It 
would be interesting to see how one could merge the 
two ends of the evolution, the perturbative with the non- 
perturbative ones. 
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